Generalized Toffoli gates using qudit catalysis 
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We present quantum networks for a n-qubit controlled gate C n ~ l {U) which use a higher dimensional (qudit) 
ancilla as a catalyser. In its simplest form the network has only n two-particle gates (qubit-qudit)- this is the 
minimum number of two-body interactions needed to couple all n + 1 subsystems (n qubits plus one ancilla). 
This class of controlled gates includes the generalised Toffoli gate C n (X) on n qubits, which plays an 
important role in several quantum algorithms and error correction. A particular example implementing this 
model is given by the dispersive limit of a generalised Jaynes-Cummings Hamiltonian of an effective spin-s 
interacting with a cavity mode. 

PACS numbers: 03.67.Lx, 03.67.Mn 



I. INTRODUCTION 

During the past decade quantum information processing 
(QIP) has changed our paradigm of computation [1]. The 
main insight was that information is physical, and hence com- 
putation, as a physical process, is limited by the laws of 
physics Quantum mechanics provides a different com- 
putational model by employing features not present in a clas- 
sical world, most notably superposition and entanglement. 
The quantum revolution opened new possibilities by devising 
novel algorifhmsj[3l 0] and tasks which either have no classi- 
cal analog IHBtzl], or provide (in some cases) an exponential 
speed-up over classically known algorithms yfl. 

A critical issue that needs to be addressed is how to con- 
struct large scale quantum circuits. As classical computing is 
designed around irreversible gates, one cannot directly trans- 
late this expertise into the quantum world. However, re- 
versible logic has a long history (8l |9[] and a central role in 
this field is played by the Toffoli gate. This gate, a controlled- 
controlled-NOT acting on three bits, flips the target bit if both 
controls are 1 . There are design tools that allow to construct 
large oracles with reversible gates. 

The Toffoli gate is also of interest in the quantum world. 
Together with the Hadamard it forms a universal set of quan- 
tum gates 1 1 Oil . Moreover, the Toffoli gate is a central build- 
ing block in phase estimation 111 ill . Shor's algorithm [3], error 
correction [ 12] and fault tolerant quantum circuits [13]. These 
properties naturally bring up the question of how to construct 
it efficiently. Since the Toffoli involves 3-body interactions, 
it does not appear naturally in physical systems, with Hamil- 
tonians usually containing 2-body interactions. One way of 
constructing the Toffoli is to decompose it into single and two 
qubit gates, the simplest decomposition in terms of CNOT 
gates require 6 such gates plus 10 single qubit gates; equiv- 
alently, one can use five two-qubit gates, but in this case these 
are general controlled-C/ instead of CNOT 11511 . One can gen- 
eralize Toffoli gates to n qubits, n — 1 controls and a tar- 
get, where now the target qubit is flipped if all the controls 
are 1. For generalized Toffoli gates the resources increase 
rapidly, requiring 0(n 2 ) two-qubit gates II 1 511 . Experimen- 
tally, 3-qubit Toffoli gates have been implemented in NMR 
systems lll2ll and ion traps ifTill . 



Recently a new approach has provided a fresh insight of 
how to construct efficiently Toff oli g ates llrjl Wh using fewer 
resources than previous designs lflll[l5tl . The key element in 
this construction Jill is the use of a larger Hilbert space by 
transforming the target qubit into a qudit. Intuitively, expand- 
ing the two dimensional space of a qubit to a n dimensional 
qudit space improves the "manoeuvrability" in the extended 
Hilbert space and results in a simplified quantum network with 
fewer gates. However, this requires that during the gate oper- 
ation the state of the system goes out of the computational 
space - effectively one of the computational qubits becomes 
a qudit. Therefore, at the end of the gate, this qudit has to 
re-enter the qubit space in order to continue the computation 
and hence needs to "forget" the non-computational degrees of 
freedom. A consequence is the potential leakage out of the 
qubit space: imperfect gates will result in leakage, especially 
if such a procedure is reiterated during a complex algorithm. 

In this article we build on this insig ht d [H but adopt 
a different approach resulting in a more efficient design. We 
completely decouple the qudit space from the computational 
space by using an ancilla. This has several important conse- 
quences. First, there is no need to expand one of the compu- 
tational qubits to a qudit and then back to a qubit. Second, 
this separation of the ancilla from the computational degrees 
of freedom allows us to measure and discard the ancilla at the 
end of the gate. Consequently we recover the computational 
Hilbert space in a pristine form - the ancilla starts and ends 
in a factorized state with respect to the qubits, hence there is 
no leakage. Here measurement plays also the role of a simple 
error correction mechanism. And finally, we have the free- 
dom to implement the qudit in a different physical system 
than the computational qubits. This means we can optimize 
other properties of the ancilla compared to the qubits (e.g., in- 
teraction strength, higher dimensional Hilbert space etc). For 
example, the computational qubits can be photons (low de- 
coherence) and the qudit ancilla can be an atom in a cavity 
(higher decoherence but stronger interaction); in this case the 
ancilla has to be coherent only for the duration of the Toffoli 
gate, not for the duration of the whole algorithm. The ability 
to measure the ancilla proves to be crucial and enables us to re- 
duce the number of entangling gates: our second design needs 
only n two-particle gates, compared to 2n — 3 in Ref . Il6l \17n . 
Note that n is the minimum number of two-body interactions 
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required to couple n + 1 subsystems (n qubits plus one an- 
cilla): less than this, at least one of them is left uncoupled 
with the rest of the system. 

We will discuss two related designs. The first does not re- 
quire measurement, but needs 2n — 1 two-body gates. The 
second design uses only n two-particle gates plus measure- 
ment of the ancilla and a feed-forward correction (single qubit 
phase shifts on the control qubits). It is important to stress that 
both networks are deterministic and the result of the measure- 
ment (in the second case) is not used for post-selection. Here 
measurement plays the same role as in teleportation: it pro- 
vides the information needed to correct the final state - the 
gate works with unit probability in principle. Finally, we give 
an example of a particular model which can be used to imple- 
ment our scheme. This corresponds to the dispersive limit of 
a generalized Jaynes-Cummings Hamiltonian of an effective 
spin-s interacting with a cavity mode. 



II. QUANTUM NETWORK FOR TOFFOLI GATE 
A. Background and notations 
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FIG. 1: An n-qubit controlled U gate C n_1 ([/); if U = X this cor- 
responds to a generalized Toffoli gate. Qubits i\ , . . . , i„ are denoted 
by thin, black lines; the ancilla (red, thick line) is a qudit with a n 
dimensional Hilbert space. The controlled-C/ gate in the middle be- 
tween a qubit and the qudit ancilla is performed only if the ancilla is 



in the |1) state, i.e., 
below the gate). 



C(U)\i n )\j) = [U"^\i n )]\j) (symbolized by 1 




I2> -X 
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Before discussing the main result of the article, let us briefly 
review a few background ideas and notations. A qudit is a n 
level quantum system; let {|0), . . . , \n — 1)} be its n dimen- 
sional (computational) basis. The action of the generalized 
Pauli operators X n and Z n in this basis is: 

X n \k) = \k + l) 

Z n \k) = u k \k) (1) 

with lj = e 2 ™/"; the addition in \k + 1} is modulo n. Thus 
Z n = diag (1, lj, lj 2 , . . . , w"" 1 ) in this basis. For n = 2, we 
denote the usual Pauli operators acting on a qubit by X and 
Z. The Fourier transform over the qudit space is defined as 

n-l 

F\k)=n- 1 / 2 J2 0jki \fi O 

3=0 

The Fourier transform maps the computational Z„-basis to the 
X n -basis (up to a relabeling); we also have FXnF^ 1 = Z n . 
For a qubit (n = 2) F is the Hadamard gate H. 

A quantum gate C n_1 (U) is defined to have n — 1 control 
qubits and performs a conditional U gate on the rt-th qubit (the 
target) if and only if all the controls are 1, i.e., .-z n _i = 1, 
with if. =0,1, The action is given by: 

C n - 1 (U)\i 1 ...i n ) = \ii...i n - 1 )U h -- i ^\i n ) (3) 

The case U — X corresponds to the generalized Toffoli gate 
acting on n qubits, namely C n ~ 1 (X); its action is to flip the 
target qubit if all the n — 1 control qubits are in the |1) state. 



FIG. 2: A measurement based C" I_1 (C7) gate. The n — 1 gates 
C(X n ) used to disentangle the ancilla from the qubits are replaced 
by a measurement of the ancilla in the Fourier basis (gate F on the 
qudit) and feed-forward. The gates P a = diag (1, ui a ) are condition- 
ally applied depending the value a of the measured ancilla; the same 
gate P a is applied to all control qubits, resulting in a homogeneous 
design. 



qudit as an ancilla. 

Consider the quantum network in Fig. Q] The controlled- 
U gate in the middle between the qudit ancilla (control) and 
i n (target) is performed only if the qudit is in the |1) state, 
i.e., \i n )\j)— >[U Sjl \i n )]\j)- With this observation, the gate 
operation becomes straightforward. For basis states, the an- 
cilla acts like a counter, adding all the values of the control 
qubits. Since the initial state of the ancilla qudit is |2), the 
middle C(U) gate is performed when 2 + i% + . . . + i n -\ = 1 
mod n, that is all control qubits are in the state The last 
n — 1 gates C(X n ) after the C(U) gate disentangle the ancilla 
from the computational qubits and return it to its initial state. 

The choice of the initial state |2) for the ancilla is deter- 
mined by the fact that the middle C(U) gate (between the 
ancilla and the nth qubit) acts only if the ancilla is in the |1) 
state. If, experimentally, we want to have the C(U) gate con- 
trolled by the \k) value, one has simply to prepare the ancilla 
in the \k + 1) state. 



C. Disentangling by measurement: a simpler design 



B. A first gate design 

Now we are ready to discuss the quantum network for a 
generalized controlled gate C" _1 (C/) using a n-dimensional 



In the previous scheme each control qubit has to interact 
twice with the qudit, first to entangle, then to disentangle it 
from the ancilla (see Fig.[T]). If the qubits are photons and the 
ancilla is an atom in a cavity, as in the photonic module lfl8l 
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[l9ll . then all the control photons have to pass twice through the 
cavity. From an experimental point of view this will require 
to store the photons in a memory or buffer and then redirect 
them back to the cavity, a situation far from optimal. Is there 
a way to simplify the network such that each qubit interacts 
only once with the ancilla? 

In Fig. |2] we present an alternative design. By measuring 
the ancilla in the Fourier basis (i.e., by first applying a Fourier 
transform F followed by a measurement in the computational 
basis) and using feed-forward, we can effectively disentangle 
the ancilla from the rest of the qubits. 

Assume we have an arbitrary input state \tpo) = 
XL- cni\i\ . . . i n )|2), where the sum is over all basis states 
i = (ii,... ,i n ) of n qubits; the ancilla is factorized and 
starts in the state |2). After applying the n — 1 controlled- 
X n gates C{X n ) on the ancilla, the state becomes \ipi) = 
J2i ctj |*i . . .i n )\ii + ■ ■ ■ + i n -i + 2} and the ancilla is now 
entangled with all the control qubits i\, . . . , i n -i- By apply- 
ing the C(U) gate between the ancilla and the target qubit i n , 
the state of the system changes to 

= J^Oilii • • • in^U 11 --^- 1 \i n )\h + . . . + i„_i + 2) 

i 

After the Fourier transform F on the ancilla, eq.©, the previ- 
ous state becomes (up to normalization): 

n-l 

Oi\ii ■ ■ ■ in-i)^ 1 --'- 1 \i n ) J2 UJ k ^ + - + ^-^\k) 
i k=0 

Suppose we now measure the ancilla and we obtain the value 
a, with < a < n — 1; the system is then projected to the 
state 

<*i\il ■ ■ • in-i)^' 1 --*"- 1 | in )^ (ll + - +l "- 1+2 V) 

i 

Hence by applying a corrective phase shift P a to all n — 1 
control qubits, with P = diag (1, uj), we obtain 

1^) = u? a Y oti\i x . . . U}^-"*- 1 \i n )\a) 

i 

which is indeed the desired state, modulo an (irrelevant) over- 
all phase uj 2a . 

It is important to stress that although the ancilla is mea- 
sured, there is no post-selection, hence the gate is determinis- 
tic and works with unit probability. Measurement here plays 
the same role as in teleportation - it provides the information 
we need to correct the state in the end. 



D. A generalized Jaynes-Cummings model 

So far the discussion was rather general without consider- 
ing a particular physical model. We now examine a possi- 
ble implementation of the abstract scheme presented above. 
The first question we ask is: What interaction Hamiltonian 
gives us the controlled unitary transformation C(X n )? Us- 
ing the identity X n = F~ 1 Z n F, the problem reduces to find 




FIG. 3: A quantum network equivalent to Fig. [2] using the identity 

X n = F 1 Z n F. 



a Hamiltonian whose evolution enacts a controlled-Z„ gate 
C{Z n ), as in Fig. El 

Assume we have a qudit (e.g., an atom in a cavity or a nu- 
clear spin in a NMR setup) interacting with a photonic qubit 
according to 

H I = X ha ji aS z (4) 

where (a) are photonic creation (annihilation) operators and 
S z = diag (s, s— 1, . . . , — s+1, — s) is the effective z-spin op- 
erator associated to the qudit, with n = 2s + 1 the dimension 
of its Hilbert space. The unitary operator induced by this (time 
independent) Hamiltonian acting for a time t on the system is 

U(t) = e - lHlt/h = exp(-i X t(JaS z ) (5) 

Let \k)\m) be a basis state of the system, with \k),k = 0, 1 a 
Fock state of the field (photon) and |m), < m < n — 1, a 
basis for the qudit. The action on this basis is 

U(t)\k)\m) = \k)[U k (t)\m)} (6) 

with U(t) = exp(—ixtS z ) a single qudit unitary. Taking 
\t = 2tt /n and recalling that ui = e 27 ™/^ we obtain 

U{2ir/ X n) = -w^diag (1, w, . . . , co^ 1 ) = -uJ 1/2 Z n (7) 

Therefore the Hamiltonian (0]) implements, up to a phase, the 
unitary we need. 

It is well-known that in the case of a qubit (n — 2) the 
Hamiltonian Hj can be obtained as the dispersive limit of the 
Jaynes-Cummings Hamiltonian Hjc = hg(a'(T- + aa + ). 
This Jaynes-Cummings Hamiltonian ll20l 12111 describes the 
interaction of an effective spin-1/2 (e.g., a 2-level atom in a 
cavity) with a quantized field. 

It is easy to show that the dispersive limit holds also for 
the general case general spin-s. Assume we have a spin-s 
(or equivalent) in a cavity interacting with a quantized field 
according to: 

H = fwa ! a + W,S z + 2hg{a ! S- + aS + ) (8) 

where u> (f2) is the resonant frequency of the cavity field 
(spin) and S± = S x ± iS y are the ladder operators for a 
spin s, with n = 2s + 1. These obey the usual commuta- 
tion relations [S' + ,S , _] = 2S Z and [S Z ,S±] = ±S±. The 
dispersive limit g/A <C 1, with A = £1 — u> the detun- 
ing, is obtained by performing the unitary transformation l22tl 
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V = exp[2g(aS+ — a^S^ )/A] and expanding to second order 
in g 

H' = VHV ] = huja f a + MIS Z + ^-{S+S- + 2a f aS z ) 
+ 0(g 3 /A 2 ) (9) 

This points towards two possible routes for physical imple- 
mentation. One possibility is the dispersive limit (g/A <C 1) 
of an actual spin s in a cavity, interacting with the field through 
the Hamiltonian of eq. ^ and with capability for S z measure- 
ment. An alternative is to use n — 1 dispersive qubits, each at 
the same frequency ft, not coupled to each other but all cou- 
pled to the field with the same strength ll23ll . With 

n— 1 n—1 
3=1 3=1 

the required commutation rules hold for the effective spin s 
(where n = 2s + 1) and the Hamiltonian of eq. ^ emerges. 
Clearly here the full (n — l)-qubit Hilbert space is bigger 
than that required. However, since the evolution preserves the 
symmetry, if the qubits are initialised in the symmetric (un- 
der qubit interchange) subspace they remain in it and act as 
an effective spin s. For the two-qubit case n = 3, the sym- 
metric states, written in terms of the individual a z eigenstates 
would be {|Tt>.^(IU) + lit)), III}}- Given the symmet- 
ric subspace, measurement projecting into the S z basis can be 
effected by measuring the individual qubits in their a z basis. 

The other core ingredient of Fig.|3]is the Fourier transform 
(FT) over the qudit. It can be implemented efficiently with at 
most n 2 — 1 pulses - any unitary U € SU (n) can be written as 

JH U = UkL'i 1 e if}kTk , with T k the SU(n) generators and 
(3k £ 1R. Proposals to implement the FT include multilevel 
atoms i25ll and harmonic oscillators l26ll . 



III. CONCLUSIONS 

A challenge for the present day quantum engineers is to 
minimize the quantum resources required to perform a given 
quantum task. In this article we discussed two quantum 
networks implementing a family of controlled gates on n 
qubits C n ^ 1 (U). This includes the generalized Toffoli gate 
C n ^ 1 (X), an important ingredient of quantum error correc- 
tion algorithms. 

Recently several authors proposed a new approach for an 
efficient implementation of generalized Toffoli gates lfl^[l7ll . 
The insight was to use a higher dimensional Hilbert space by 
allowing one of the qubits (the target) to temporarily become 
a qudit during the gate operation. The end result is a reduction 



in the number of two-body quantum gates required at the ex- 
pense of going out and in of the computational space during 
gate execution. 

Let us now review the main features of our construction. 
Our design strategy was to decouple the higher dimensional 
space (the qudit) from the computational qubits by choosing 
the qudit to be an ancilla. This has several consequences: first, 
it prevents leakage from the computational space, as the an- 
cilla can be discarded at the end of the gate. Second, the qubits 
and the ancilla can be implemented in different physical sys- 
tems (e.g., photons vs. an atom in a cavity). And third, as the 
ancilla does not have any computational degrees of freedom, 
it can be measured without destroying the coherence of the 
qubits. This is the case of the gate design in Fig. [2] where 
measuring the ancilla in the Fourier basis disentangles it from 
the computational qubits and simplifies the design: each qubit 
interacts minimally with the ancilla - only once. Another ad- 
vantage of the present scheme is that all control qubits interact 
homogeneously, i.e., in the same way, with the qudit; more- 
over, there is no need of extra single qudit gates on the ancilla 
between the C(X n ) gates (in contrast, in Ref. lfT6l [l7ll extra 
single-qudit gates are required). This last property is espe- 
cially useful for quantum computation schemes in which the 
ancilla is symmetrically coupled to all control qubits. One 
particular example is a NMR approach where a central qubit 
is coupled identically (due to symmetry) to a number of satel- 
lite qubits, as in Ref. [27]. In this case we can enact the n—1 
C(X n ) gates between the control qubits and the ancilla with 
a single NMR pulse, provided the target qubit is distinct from 
the rest, hence has a different coupling. 

The controlled C(Z n ) unitary we use in our scheme can 
be implemented with a Hamiltonian which generalizes the 
dispersive limit of Jaynes-Cummings interaction of a qudit 
coupled to a photonic field. One physical implementation of 
this system could be an actual matter qudit (a spin s, with 
n = 2s + 1) coupled to a cavity field. Another possibility is 
to realise an effective qudit from n—1 qubits each coupled to 
the field. The latter approach could be feasible with supercon- 
ducting qubits each positioned at an antinode of a microwave 
cavity field [28]. For this example the required measurement 
projecting into the S z basis might be effected through disper- 
sive measurement of the individual qubits |29[ l30ll at their op- 
erating points, or possibly through gate voltage bias change 
followed by charge measurement. 
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